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ABSTRACT 

Spider is a long-duration, balloon-borne polarimeter designed to measure large scale Cosmic Microwave 
Background (CMB) polarization with very high sensitivity and control of systematics. The instrument 
will map over half the sky with degree angular resolution in I, Q and U Stokes parameters, in four 
frequency bands from 96 to 275 GHz. Spider^s ultimate goal is to detect the primordial gravity wave 
signal imprinted on the CMB B-mode polarization. One of the challenges in achieving this goal is the 
minimization of the contamination of B-modes by systematic effects. This paper explores a number of 
instrument systematics and observing strategies in order to optimize B-mode sensitivity. This is done by 
injecting realistic-amplitude, time-varying systematics in a set of simulated time-streams. Tests of the 
impact of detector noise characteristics, pointing jitter, payload pendulations, polarization angle offsets, 
beam systematics and receiver gain drifts are shown. Spider^s default observing strategy is to spin 
continuously in azimuth, with polarization modulation achieved by either a rapidly spinning half-wave 
plate or a rapidly spinning gondola and a slowly stepped half-wave plate. Although the latter is more 
susceptible to systematics, results shown here indicate that either mode of operation can be used by 
Spider. 

Subject headings: cosmic microwave background, polarization experiments, B-modes, gravity waves, 
analytical methods 



1. INTRODUCTION 

In the past decade, a wealth of data have pointed to 
a "standard model" of the Universe, composed of ~ 5% 
ordinary matter, ~ 22% dark matter and ~ 73% dark en- 
ergy in a flat geometry. The flatness of the Universe, the 
near isotropy of the CMB, and the nearly-scale-invariant 
nature of the primordial scalar perturbations from which 
structure grew support the existence of an early acceler- 
ating phase dubbed "inflation". A necessary by-product 
of inflation is tensor perturbations from quantum fluctu- 
ations in gravity waves. A detection of this Cosmological 
Gravity- Wave Background (CGB) would give strong evi- 
dence of an inflationary period and determine its energy 
scale, while a powerful upper limit would point to more 
radical inflationary scenarios, e.g., involving string theory, 
or some non-inflationary explanation of the observations. 

The CGB imprints a unique signal in the curl-like, or 
B-mode, component of the polarization of the CMB; de- 
tection of a B-mode signal can be used to infer the presence 



of a CGB at the time of decoupling. Direct detection of 
the gravity waves is many decades off; an advanced Big 
Bang Observer success or to LISA h as been suggeste d as 
a wa y to achieve this (jPhinnev et a l. 2004; Harry e t al.l 
l2006t ). Thus a measurement of the primordial B-modes is 
the only feasible near-term way to detect the CGB and 
have a new wind ow to the physics of the early Universe 
(jBock et al.ll20Q6[ ). 

A CGB with a potentially measurable amplitude is 
a by-product of the simplest models of single field in- 
flation which can reproduce the scalar sp ectral tilt ob- 
serve d in current coinbined CMB data (jSpergel et al.l 
120071: iMacTavish et al.]|2006D . Examples ar e chaotic in- 
flation froin pow er law inflaton potentials (|Lindd 1X9831 : 
iLinde et al.l |2005[) or natural inflation from cosine infla- 
ton pote ntials involving a ngular (axionic) degrees of free- 
dom (A dams et al.lll993[) . The amplitude is usually pa- 
rameterized in terms of the ratio of the tensor power spec- 
trum to the scalar power spectrum, r = 'Pt{kp)/Vs{kp), 
evaluated at a comoving wavenumber pivot kp, typically 
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taken to be 0.002 Mpc^ . Chaotic inflation predicts 
r w 0.13 for a cjP potential and r « 0.26 for a potential, 
and natural inflation predicts r « 0.02 — 0.05. 

The potential energy V driving inflation is related to 
r by y « (10^^ Gev)'*r/0.1. Low energy inflation mod- 
els have low or negligible amplitudes for the CGB. To get 
the observed scalar slope and yet small r requires special 
tuning of the potential. These are often more complicated, 
multiple field models, e.g., (Linde 1994), or string-inspired 
brane or moduli models (Kallosh 2003). Given the collec- 
tion of models it is difficult to predict a precise range for 
the expected tensor level and the prior probability for r 
should be considered as wide open. 

Recent CMB data have reached the sensitivity level 
required to constrain the amplitude of and possibly 
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nificant complication of the 



measurements is that the E-mode amplitude is an order 
of magnitude lower than the total intensity. In addition, 
galactic foregrounds such as synchrotron and dust are ex- 
pecte d to be significant at these amplitudes (|Kogut et al.l 
l2007t ). Furthermore the polarization properties of fore- 
grounds are largely unknown. Constraining B-modes 
presents an even greater challenge as it is a near certainty 
that polarized foregrounds will dominate the signal. 

The next generation of CMB experiments will benefit 
from a revolution in detector fabricat ion in the form of 
arrays of antenn a-coupled bolometers (iGoldin et al.ll206^ 
iKuo et aTll2006D . The antenna-coupled design is entirely 
photo-lithographically fabricated, greatly simplifying de- 
tector production. In addition, the densely populated an- 
tennas allow a very efficient use of the focal plane area. 

Spider will make use of this technological advance, in 
the form of 2624 polarization sensitive detectors observing 
in four frequency bands from 96 to 275 GHz. A multi- 
frequency observing strategy is a necessary requirement 
to allow for a subtraction of the foreground signal. Spi- 
der will observe over a large fraction of the sky at degree 
scale resolution producing high signal-to-noise polarization 
maps of the foregrounds at each frequency. 

Extraordinarily precise understanding of systematic ef- 
fects within the telescope will be required to measure the 
tiny B-mode signal. This paper presents a detailed investi- 
gation of experimental effects which may impact Spider^s 
measurement of B-modes. The strategy is to simulate a 
Spider flight time-stream injecting systematic effects in 
the time domain. The aim is to determine the level of B- 
mode contamination at subsequent stages of the analysis. 
With these results one can set stringent requirements on 
experimental design criteria, in addition to optimizing the 
telescope's observing strategy. 

The outline of this paper is as follows. Section[2]gives an 
overview of the instrument, flight and observing strategy. 
Section [3] describes the details of the simulation method- 
ology. Results for several systematic effects are presented 
in Section [J] Section [5] concludes with a summary and 
discussion of the results. 

2. THE INSTRUMENT 



An initial descriptio n of Spider can be found 
in iMontrov et al.l (|2006b| ). Since that publication some of 
the telescope features have been changed in order to sim- 
plify design and to further optimize the instrument. An 
overview of Spider instrumentation is given here. 

Spider is a balloon-borne polarimeter designed to mea- 
sure the polarization of the CMB at large angular scales. 
Beyond the ability to map large scales, a clear advantage of 
a balloon platform is the increase in raw sensitivity, espe- 
cially in the higher frequency channels, achievable above 
the Earth's atmosphere. For the long-duration balloon 
(LDB) flight Spider will launch from Australia, with a 
^--^25 day, around the world, constant latitude, trajectory. 
The first test flight is scheduled for fall 2009, a turn-around 
flight from Alice Springs, Australia, of ^ 48 hours dura- 
tion. 

A schematic of the Spider payload is depicted in the 
left panel of Figure [1] The Spider gondola will spin in 
azimuth at a fixed elevation, observing only when the sun 
is a few degrees below the horizon H A constant latitude 
25-day flight, launching from Australia (with the optical 
axis tilted at 41 degrees from the Zenith), yields a sky 
coverage of ~ 60%. 

Azimuthal attitude control is provided by a reaction 
wheel below the payload and by a torque motor in the 
pivot located above the gondola. Spider will employ a 
number of sensors to obtain both short and long time 
scale pointing solutions. These include: two star cam- 
eras, 3 gyroscopes, a GPS and a three axis magnetometer. 
The pointing syst em is based on proven BOOMER ANG 
(jMasi et al.ll2006l ) and BLAST (jPascale et al.ll2007D tech- 
niques. The pair of star cameras are mounted above the 
cryostat on a rotating platform, which will allow them to 
remain fixed on the sky, providing pointing reconstruction 
accurate to ~ 6". Solar arrays pointing toward the sun 
during daytime operation will recharge the batteries sup- 
plying payload power. 

The instrument consists of 6 monochromatic telescopes 
operating from 96 GHz up to 275 GHz. All six telescope 
inserts are housed in a single LHe cryostat which provides 
>30 days of cooling power at 4K (for the optics) and at 
1.5K (for the sub-K cooler). The detectors are further 
cooled to 300mK using simple '^He closed cycle sorption 
fridges, one per insert, which are cycled each day when 
the sun prevents observations. Specifications for each of 
the 6 telescopes, including observing bands and detector 
sensitivities are given in Table [T] 

Spider uses antenna-cou pled bolometer arrays cooled to 
300 mK (|Kuo et al.ll200il . Figure [2] shows an image of a 
prototype detector and the measured beam response of a 
single dual-polarization antenna. The antenna arrays are 
intrinsically polarization sensitive, with highly symmetri- 
cal beams on the sky and low sidelobes. Each spatial pixel 
consists of phased array of 288 slot dipole antennas, with 
a radiation pattern defined by the coherent interference 
of the antennas elements. Each of the feed antennas pro- 
vides an edge taper of roughly — 13 ± 1 dB on the primary 
aperture. A single spatial pixel has orthogonally polar- 
ized antennas. The optical power incident on an antenna 
is transmitted to a bolometer and detected with a super- 
conducting transition-edge sensor (TES) immediately ad- 



An additional daytime (anti-sun) scanning mode may be implemented but is not discussed here. 
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Fig. 1. — Left: The Spider payload. Six independent monochromatic telescopes are housed in a single long hold time cryostat. 
Each telescope is fully baffled from radiation from the ground. Power is supplied by solar arrays. The baseline observing strategy 
is to spin the payload in azimuth at fixed elevation. Spider is designed to obtain maximum sky coverage during a 20-30 day, 
mid-latitude, around-the-world flight. Right: Spider optical train. The telescope yields a flat and telecentric focal plane. The 
apodized Lyot-stop, which is fixed with regard to the instrument, is maintained at 4 Kelvin. All dimensions are in millimeters. 



jacent to the spatial pixel. 

The TES sensors will be read out using superconducting 
quantum interference device (SQUI D) current amplifiers 
with time-domain multiplexing (.Chervenak et al.l 199^ 



de Korte etall 120031 : iReintsema et al.l 120031 : llrwin et al.l 

Ambient temperature multi-channel electronics 



(MCE) (jBattisteUi et al. 2007), initially developed for 
2II0 



shown in the 
on the Robin- 

The monochro- 
two AR-coated 



SCUBA2 (|Holland et al.ll2006( ). will work in concert with 
the time-domain multiplexers. 

The optical design for the inserts 
right panel of Figure [Tj is based 
son /BICEP (jKentins et al^T2003^ optics, 
matic, telecentric refractor comprises 
polyethylene lenses and is cooled to 4K in order to reduce 
the instrumental background to negligible levels. The pri- 
mary optic is 302 mm in diameter and the clear aperture 
of the Lyot-stop is 264 mm which produces a 45' beam at 
145 GHz. 

A cryogenic half-wave plate is located in front of the 
Lyot-stop of each telescope. Rotating the half-wave plate 
aids in polarization modulation, making for cleaner mea- 
surements of the Stokes Q and U, mitigating the need to 
difference detector time-streams. This is essential to the 
reduction of the requirements for characterization of indi- 
vidual detectors and ultimately a reduction of experimen- 
tal systematics. The half-wave plate consists of a single 
birefringent sapphire plate coated with a single layer of 
Herasil quartz on each side. 

Initially polarization modulation was to be achie ved via 
a con tinuously spinning half-wave plate (Montro v et al.l 
120061)1 ) . This work examines the viability of a fast spinning 



gondola modulating the incoming signal with the half- wave 
plate stepping 22.5 degrees per day. Section I4d1 illustrates 
that either of these modes of operation can be used for 
Spider. The latter mode, stepping the half-wave plate, is 
easier to design mechanically and more robust to operate, 
and is therefore preferred. 

3. SIMULATION METHODOLOGY 

The simulation pipeline is based la r gely o n the anal- 
ysis pipeline described in I Jones et"al] (|2007f) which was 
developed for the analysis of the data obtained from ob- 
servations made with the BOOMERANG telescope dur- 
ing th e 2003 LDB Antarctic fligh t (Montrov et al. 2006at 
iPiacentini et al.ll2006l : I Jones et aL|[2006l : lMasi et al.ll2006[ ) 



A schematic outlining the components of the simulation 
pipeline and the various inputs and outputs is given in 
Figure [31 

The flight simulator generates time-ordered pointing 
data in the form of right ascension, declination and po- 
larization angle for each detector. Data are simulated for 
16 detectors, arranged in evenly spaced pairs in a single 
column which extends the full height of the focal plane. 
Detectors in a pair are oriented to be sensitive to orthog- 
onal polarizations. Since signal-only simulations are used 
in this work this is sufficient to test most of the system- 
atic effects and observing strategies considered here. The 
small number of time-streams also reduces signiflcantly the 
data storage and computation requirements which itself 
will present a unique challenge for the actual analysis. 

It is assumed that the telescope is flxed in elevation 
41 degrees from the zenith and that the payload is mov- 



4 



MacTavish et al. 



Obs. Band (GHz) 


96 


96 


145 


145 


225 


275 


Orientation 


Q 


U 


Q 


U 


Q 


U 


Bandwidth (GHz) 


24 


24 


35 


35 


54 


66 


Number of Detectors 


288 


288 


512 


512 


512 


512 


NET ifiKy/s) 


100 


100 


100 


100 


204 


351 


Beam FWHM (arcmin) 


58 


58 


40 


40 


26 


21 



Table 1 

Spider Channel Specifications. Instrument orientation, observing bands, detector counts, sensitivities and beams. A total of 2624 
detectors is distributed between six telescopes, with two operating at 96 GHz and two operating at 145 GHz. 




Fig. 2. — Left: A single pixel of a 145 GHz antenna-coupled bolometer, comprising a 288 element phased array of dual-polarization 
slot antennas coupled to a matched load by a superconducting microstrip network. Microstrip filters, which determine the spectral 
response, and TES detectors, which measure the power dissipated in the load are visible at bottom. Right: The measured beam 
pattern of the dual-polarization antenna. The upper limit on differential beam ellipticity of is 1%, limited by the testbed. The 
polarization efficiency is greater than 98%. It is important to note that the beam pattern here is the feed beam pattern. The 
beam on the sky is infiuenced by the telescope. While the Spider telescope edge taper is modest, the beam on the sky will be 
more symmetric than the feed pattern shown here. In particular, the visibly large and asymmetric lobes above will not propagate 
to the sky. 



ing in longitude (beginning at 128.5 east) at a speed of 
3.76 X 10"'' degrees per second at constant latitude (25.5 
south). Data are simulated for four days of operation, 
assuming a mid-November launch. Four days operation 
allows for one complete observing cycle for the stepped 
half-wave plate operating mode, after which the cycle is 
repeated. This is also the minimum required to ensure 
sufficient coverage for polarization reconstruction of the 
entire observed area. 

Full sky intensity and polarization maps are simulated 
and smoothed with the Spider beam using th e synfast pro- 
gram which is part of the HEALPix software (|G6rski et al.l 
l2005f ). In order to ensure that signal variation within a 
pixel is negligibly small the full sky maps are pixelized at 
a resolution which corresponds to a pixel size of ~ 3.4'. 

Full sky maps are then converted into time-ordered data 
(TOD) using pointing information from the flight simula- 
tor. Thus, the time-stream generator constructs di for 
each detector from the equation 

di = G[Ipix + ((3pixCOs(2V'i) + C/p„sin(2?/'i))]- (1) 
z p 

Here p parametrizes the polarization efficiency, '0 is the 
final projection of the orientation of a detector on the sky 
and G is the detector gain or responsivity. 



All time-streams are high-pass filtered at 10 niHz during 
the map-making stage. This is done to test the impact of 
the filtering that is required in the case of real data which is 
effected by long-timescale systematics. Particular system- 
atics of concern are knowledge of system transfer functions 
(or equivalently knowledge of the gains) and knowledge of 
the noise amplitude/statistics on long-timescales. 

During the time-stream generation the (stepped or spin- 
ning) half-wave plate polarization angle is added to the 
intrinsic polarization angle of the individual detectors. In 
addition polarization angle systematics, beam offsets and 
gain drift are also applied during time-stream generation. 
Additional pointing jitter and pendulation systematics are 
added to the pointing time-streams during flight simula- 
tion. For the case of multiple beam distortions, multiple 
pointing time-streams are produced and the full-sky map 
is observed by each beam. 

Finally, Spider maps are constructed in terms of the 
observed Stokes parameters, 1°^^ , Q°^^, and U°^^ with an 
iterative map maker that uses an adaptation of the Ja- 
cobi method (described in detail in I Jones et al.l (|2007t )) to 
recover the input signal. To reduce computation time out- 
put maps are pixelized at a resolution which corresponds 
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to a pixel size of ^ 13.7' (about 1/4 of the beam size). 
Although the simulations are for pure signal, the map 
maker algorithm performs an inverse noise filtering of the 
time-streams. This filtering would be included in noisy 
time-streams to reduce the strongest effects of 1// noise 
which can significantly reduce map-making efRciency. In 
this case it is also included to make the simulations used 
here accurate representations of the full pipeline. 

3.1. Residual measure 

The aim is to quantify the contamination of the B-mode 
angular power spectrum (BB) from systematics which in- 
duce either I^ Q,UorQ ^ U mixing. To assess the 
impact of the various systematics on BB the following pro- 
cedure is implemented: 

• Generate full-sky input /. Q. and U maps with 
Cf^ =0. 

• 'Observe' the maps with simulation pipeline includ- 
ing a chosen systematic (but no noise) giving signal 
only Q°^<', and U°^\ 

• Take the difference between input and output maps 
over the survey area 



Jl cs 


yobs 


-I, 






= Q°^' 


-Q, 


(2) 




^obs 


- u. 





• Apply a mask with pixel weighting determined by 
the number of observations per pixel to the residual 
maps. 

• Spherical harmonic transform the weighted maps to 
obtain pseudo-C^ spectra of the BB residual. 

• Compute residual measure i?f ^ defined below. 

The pixel weighting is applied in order to reduce the ef- 
fect of badly sampled pixels at the edge of the map. These 
would bias heavily the raw pseudo-Cf computed from the 
residual maps. The mask that is applied to all of the sim- 
ulations for which the gondola spin rate is 36 degrees per 
second (dps) is shown in Figure [H The dark band rep- 
resents a galactic cut at ±10° in galactic latitude. The 
white region represents the portion of the sky that cannot 
be observed because of the sun. With these regions flagged 
the fraction of the sky covered for this observing strategy 
is ~ 60%. Each pixel value in the mask is the number of 
observations in the pixel divided by the value in the pixel 
with the maximum hits. 

The spectra obtained from the met hod above are raw 
cut-sky, or pseudo-C^, power spectra (jHivon et al.l 12002} ). 
Since no B-mode power is present in the original full-sky 
simulation any B-mode power in the final maps will be 
due to the mixing of modes from either systematics or cut- 
sky effects which mix E and B-modes |Lewis et al.ll2002t 
iBunn et"al]|2003D . In the limit of full-sky a simple model 
for the spherical harmonic coefficients of the residual maps 
is given by 

afr - [Ff-n + [Frn aj^. (3) 



The terms on the right hand side of Eq.[3]represent leakage 
of total intensity and E-mode into B due to time domain 
effects induced by the pipeline. They include any loss of 
modes due to time domain filtering. The total transfer is 
described as an isotropic coefficient Fg . There is wo E ^ B 
mixing due to any cuts in this case. 

A significant assumption introduced above is that time 
domain effects result in isotropized transfer of mode power 
in the map do main, hence no m -dependence of the trans- 
fer functions (jHivon et al.ll200^ . The validity of the as- 
sumption depends on the observing strategy adopted but 
should be appropriate for polarization experiments where 
cross-linking is maximized. 

When including cut-sky effects the pseudo-C^ residual 
B-mode spectrum can then be approximated as 

2 



c, 



BBres (noBB) _ 



1/2 



- 1 



EE 



' M, 



r2 p(2)^TT 



7^ n2 p(3)r'TE 



(4) 



Here the Mu' is the geometric total intensity kernel, K„, 
is the geometric leakage kernel (Szapudi et al. 200H) cou- 
pling E ^ B, and xK/,^, is the geometric kernel for the 

cross-correlation spectrum. The transfer functions F^}'"^'^^ 
represent a combination of individual transfer effects 



- 

^ = ( \Fn 



The quantity, C, 



C 



BBres (noBB) 



1/2 



.T^Bl 1/2 



(5) 
(6) 

(7) 



in dH) is divided by 
^BB(noEE) ^-^^ gg signal obtained from reconstructed 
Spider Q/U maps with only an input BB signal (for r = 
0.1) and no input EE 

^BB{noEE) x,, j,2 j^BBr^BB 

- 2^+Ku'B^,I<^, Cg, . 



(8) 



Dividing ([1]) by ([5]) gives the fraction of the raw BB power 
that is coming from a transfer effect caused by systematics 
and not from an input BB spectrum. The same mask, with 
pixel weighting determined by the number of observations 
per pixel, is applied to all maps (residual and Q/U) which 
have the same gondola spin rate. The final BB residual 
measure is then defined as 

^BBres (noBB) 

pSB _ r<BB fa\ 

- p^BB (noEE) ' '^^) 

This multiplies the fractional residual by the input BB 
spectrum for r = 0.1 giving the residual in terms of an 
equivalent BB signal. The residual measure defined above 
is not designed to give a complete picture of how well the 
original BB signal can be reconstructed from the observa- 
tions. A complete treatment would require a full un-biased 
power spectrum estimation method, which is beyond the 
scope of this work. Instead (O isolates the impact of the 
systematics under study on the observed signal by min- 
imizing the impact of the E B mixing from cut-sky 
effects. 



6 



MacTavish et al. 



Sky ^ £ 

Model 




Systematics 

->Noise Filter 
->Pol. angle offsets 
->Multi-beam TODs 
->Beam offsets 
->Gair, drift 



Sky Simulator 
(Synfast) 



■ 


Full Sky 






IQU 




■ 


Maps 


■ 



I 



Gain ) (^FU^^ 



Timestream Generator 



Systematics 
->Pointing jitter 
->Pendulations 



Flight Simulator 



7^ 




Spider 
Observed 
IQU 
Maps 



Systematics 
->Mulit-beam pointing 



Fig. 3. — A schematic representation of the simulation pipeline. During the time-stream generation the (stepped or spinning) 
half-wave plate polarization angle is added to the intrinsic polarization angle of the individual detectors. In addition polarization 
angle systematics, beam offsets and gain drift are also applied during time-stream generation. Pointing jitter and pendulation 
systematics are added to the pointing time-streams during flight simulation. For the case of multiple beam distortions, multiple 
pointing time-streams are produced and the full-sky map is observed by each beam/pointing. 
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Fig. 4. — Left: The Spider mask for 36 dps simulations. The dark band represents a galactic cut at ±10° in galactic latitude and 
the white region represents sun flagging. With these regions flagged the fraction of the sky covered for this observing strategy is 
~ 60%. Pixel values are the number of observations in the pixel divided by the maximum hits value. The most obvious features are 
constant declination lines where scan circles on the sky for each detector overlap and the coverage is deepest. The pixel weighting 
is applied in order to reduce the effect of badly sampled pixels at the edge of the ma p. Top Right: Spider coverage projected into 
equatorial coordinates. Bottom Right: The IRAS 100 pm map (jSchlegel et al.|[l998l ) of Galactic dust is shown for comparison. 



Note that for all of the input maps the same initial seed 
value is used to generate the full CMB sky, i.e. the sample 
scatter is the same for all simulations. 



4. SIMULATION RESULTS 

The presentation of results begins by illustrating the 
base residuals for two basic modes of half-wave plate 
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operation-stepped and continuously rotating. For the re- 
maining subsections the B-mode residuals from experi- 
mental systematic effects for the stepped half-wave plate 
case are examined. All simulations are for signal only 
(with no noise) but time-streams are inverse noise filtered 
during the map making phase. Aside from Section 14.21 
which explores two knee frequency values, the 1/f knee 
for the noise filter is 100 mHz for all simulations. In all 
plots the case labeled nominal is a 36 dps gondola spin 
rate, with the half-wave plate stepping 22.5° once per day, 
with 10 iterations (sufficient to recover the residual levels 
of the continuously-rotating half-wav e plate case) of t he 
map-maker, a Jacobi iterative solver (jjones et al.ll2007fl . 

4.1. Polarization Modulation 

Since Spider'' s default observing strategy is to spin con- 
tinuously in azimuth, two modes of polarization modula- 
tion are explored. For the first case the half-wave plate 
spins continuously at 10 Hz while the gondola rotates at 6 
dps. For the second case the half- wave plate is stepped by 
22.5 degrees per day, while the gondola rotates at 36dps or 
more. Therefore in the first instance the half-wave plate 
is modulating the incoming polarization signal and in the 
second case the gondola itself is used to modulate the sig- 
nal. 

Modulation by the gondola spin has a number of de- 
sign advantages over the inclusion, in the optical train, 
of a continuously rotating half-wave plate. The half-wave 
plate adds a degree of complexity in the design with a 
subsequent impact on the robustness of the instrument. 
In addition it is a potential source of a number of sys- 
tematic effects for example microphonics, thermalization 
effects, magnetic pickup and higher power dissipation at 
4K. It is therefore preferable (and nearly equivalent as will 
be shown) to step the half-wave plate once per day, in or- 
der to increase Q and U redundancy in a single pixel, while 
rapidly spinning the gondola in order to move the signal 
above the detector 1// knee frequency. 

Figure [5] shows Q residual maps for the two modulation 
modes. Maps of the U residuals are not shown here but 
are of similar amplitude. The top panel shows the resid- 
uals for the first case (continuous half-wave plate rotation 
at lOHz). For this case a "naive", or zero-iteration, map is 
shown. The naive map is equivalent to a simple (pixel-hit 
weighted) binning of the time-stream into pixels. Itera- 
tions of the map-making step are not required in this case 
since the signal is modulated to frequencies higher than the 
expected 1// knee and the loss of modes at low frequen- 
cies has virtually no impact on the final maps. This is one 
of the benefits of a design which includes a continuously 
rotating half-wave plate. 

The middle panel of Figure [5] shows the naive map for 
the stepped half-wave plate mode with the gondola spin- 
ning at 36 dps. In this case significant striping is present 
due to the loss of low frequency modes. These translate to 
large scale modes along the individual scans and result in 
the striping obvious in the maps. Iterating the map-maker 
in this case reduces the effect of striping as the large scale 
modes are recovered. After 10 iterations the striping is 
significantly reduced as shown in the bottom panel of Fig- 

A full explorati o n of faster, sub-optimal map-making algorithms is left for future work. In particular do-striping algorithms (see 
e.g. lAshdown eTall ^MA\ ) may provide a much faster alternative although it is still not clear that these can be applied to a Spider ob- 
serving strategy and polarization sensitivity requirements. 



ure [51 One possible way to reduce the computational load 
of a map-making stage with many iterations is to spin the 
gondola faster to modulate the signal into higher frequen- 
cies. 

The power spectra for the BB residual measure © are 
shown in Figure [Sj The optimal solution is the continuous 
half-wave plate modulation scheme. This yields the lowest 
residual compared to an r = 0.1 fiducial BB model. In the 
stepped the 36 dps, non-iterated case the residuals have 
the same amplitude as the model on the largest scales. 
The residuals are reduced to < 1% levels for multipoles 
I < 100 when the map-maker is iterated 10 times. The 
stepped 70 dps spin case with 10 iterations yields even 
smaller residuals at the largest scales. 

Given the design and implementation advantages, the 
simple stepped half-wave plate system appears a feasible 
choice for the Spider scan strategy, albeit with significant 
additional computational costs Ej- The remainder of this 
Section is restricted to the stepped half-wave plate case. In 
particular the focus will be to probe whether any other sys- 
tematic effects invalidate this choice of modulation scheme. 

4.2. Noise 

To examine the impact of different 1/f profiles on the 
stepped mode residuals a number of different cases were 
run. 

• 36 dps gondola spin rate with 100 niHz detector 1/f 
knee. 

• 36 dps gondola spin rate with 500 mHz 1/f knee. 

• 110 dps gondola spin rate with 500 mHz 1/f knee. 

The time-stream high-pass filter cut-off is kept at 10 
mHz in all cases. A comparison of BB signal residuals 
varying the detector knee frequency is shown in Figure [71 
Although simulations are signal only the time-streams are 
inverse noise filtered, as would be done for the real data. 
This reveals the impact of the detector noise characteris- 
tics in terms of the degradation of the polarization signal 
on the largest scales. The effect of a 500 mHz knee is 
clearly seen on the largest angular scales. Even for 10 it- 
erations of the map-maker the residuals are close to the 
20% level for this case. Increasing the spin rate to 110 dps 
reduces the impact of the higher 1/f knee and approaches 
the nominal 36 dps, 100 mHz knee case. A 500 mHz knee 
frequency for the detectors and readout electronics is pes- 
simistic, but would not be catastrophic since polarization 
modulation can still be achieved by the faster spinning 
gondola. Spider^s high frequency response is limited by 
the noise and response time of the detectors the combina- 
tion of which sets the maximum gondola spin rate. With 
5 ms detectors. Spider can spin up to 110 dps before being 
affected by the detector noise and time constants. Thus 
for the stepped half-wave plate case, the limit for the 1/f 
knee frequency is ^500 mHz (for r = 0.1). 



4.3. Pointing Systematics 
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Sdps Gondola Spin Rate, Fasf Spinning Waveplate flOHg}, Malve Map 




Fig. 5. — Maps of the residuals in the Q Stokes parameters for the two polarization modulation strategies. The top panel shows 
the residuals for the continuous half- wave plate rotation case. For this case a "naive" , or zero-iteration, map is shown. The middle 
panel shows the naive map for the stepped half-wave plate mode with the gondola spinning at 36 dps. In this case significant 
striping is present due to the loss of low frequency modes. This is caused by the inverse noise filtering of the time-streams during 
the map- making phase which uses a noise kernel with a realistic 100 mHz 1/f knee. For the stepped case the polarization mod- 
ulation is not sufficient. However, iterated map-making reduces the impact of the striping as shown in the bottom panel and 10 
iterations of the map-maker are sufficient to recover most of the lost modes. 



One of the more challenging aspects of balloon-borne 
telescope observations is pointing reconstruction, consti- 
tuted as one of the limiting systematics in the inter- 
pretation of CMB data obtained from a balloon-borne 
telescopes. For Spider the poi nting systein is based 
on largel y on BOO MERANG ()Masi et al.l |2006[ ) and 
BLAST (jPascale et al. 2007) instrument pointing. The 
former achieved ~ 1.5 arcminute pointing and the latter 
reconstructed pointing to a few arcseconds. The require- 
ment for Spider will be much less demanding given the 
low resolution. Nonetheless the design goal for Spider is 
sub-arcminute pointing reconstruction and establishing a 
precise requirement is still important given the particular 
sensitivity of polarization measurements to offsets in the 
pointing. 

Two main pointing offsets are explored in this work. 
The first is a random pointing jitter added to the original 
pointing solution. The offsets are added to the right as- 
cension (RA) and declination (DEC) value of each sample 
in the solution. The offsets are constant for 6 seconds, 
after which new random values are drawn. The 6 second 



bandwidth regime is chosen since this will mimic point- 
ing systematics which occur within the time-scale of one 
gondola spin 0- A nominal run with 1 arcminute RMS 
and a worst case scenario with 10 arcminute RMS for the 
instrument jitter are considered. To translate from instru- 
ment jitter to true jitter on the sky a factor of cos(DEC) 
is applied to the RA offsets. 

The second systematic consists of a sinusoidal oscilla- 
tion with an amplitude of 6 arcminutes and a 20 minute 
period. This effect simulates the pendulation of the gon- 
dola. The pointing offsets examined are typical of in-flight 
conditions-albeit the 10 arcminute jitter is extremely pes- 
simistic. 

The results for the residual measure for the three cases 
are shown in Figure [51 The pendulation case also includes 
a long-timescale 1 arcminute RMS jitter. All cases except 
the large 10 arcminute RMS jitter result in negligible con- 
tributions to the residual measure at ^ < 100 compared 
to the nominal stepped half-wave plate/36 dps spin mode 
without any systematic. 

An additional pointing systematic affecting polarization 



White-noise pointing jitter, of similar amplitude, that varied from sample-to-sample was also explored and was found to have an insignificant 
effect. 
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Fig. 6. — Comparison of the BB residual from continuously spinning and stepped half-wave plate polarization modulation schemes. 
For the first case the half-wave plate spins continuously at 10 Hz while the gondola rotates at 6 dps. For the second case the 
half-wave plate is stepped by 22.5 degrees per day, while the gondola rotates at 36 dps or more. Operation in stepped half-wave 
plate mode, with the modulation provided by the spinning gondola, requires iterated map-making as the signal is not modulated 
as far from the 1// knee as in the continuously- rotating half- wave plate mode. 10 iterations of the map-maker is sufficient to 
recover the residual levels of the continuously-rotating half-wave plate case. Fewer iterations may be required if the gondola spin 
rate is even higher (70 dps case). The input BB spectrum for r — 0.1 is plotted for comparison. 
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Fig. 7. — Impact of a higher 1// detector knee on the BB residuals. Simulations are for signal only (with no noise) but time- 
streams are inverse noise filtered in the map-maker, giving a realistic estimate of modes that would be lost from 1/f effects in the 
real data. With the higher knee frequency of 500 mHz the gondola spin rate must be increased to reduce the residuals to <2% 
levels over the target range of multipoles I < 100. The maximum gondola spin rate of ~ 110 dps is limited by the time response 
and noise characteristics of the Spider detectors. The input BB spectrum for r = 0.1 is plotted for comparison. 
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Fig. 8. — Three cases of pointing jitter. In the first two offsets are added to the RA and DEC with new offset values (with 
RMS amplitudes 1 arcminute or 10 arcminutes) every 6 seconds. In the third case, a sinusoidal oscillation is implemented with 
an amplitude of 6 arcminutes and a 20 minute period. This effect simulates the pendulation of the gondola. In addition for the 
latter case, a long time scale (day period) 1 arcminute RMS jitter is added. Both types of pointing error, reconstruction error and 
in-flight pendulations, have negligible effects; a less than 1% effect even for the large 10 arcminute jitter. 
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Fig. 9. — Residuals from absolute and relative offset of detector polarization angles. In the first two cases random 0.5 and 1 
degree RMS ^j) errors are added to each detector. Residuals for these relative offsets are negligible. In the remaining cases the 
same offset (0.25 degree, 0.5 degree and 1 degree) is applied to all channels, simulating an overall calibration error in the half-wave 
plate angle. For a 1 degree absolute offset the effect is as high as 20%. The 0.5 degree offsets are required having only a few 
percent residual effect. 



measurements is the requirement to reconstruct the angle 
if) (in equation [T|) of the detector polarization relative to 
the fixed, local Q and U frame of reference on the sky. The 
systematic can arise in two distinct ways. The first is a 
relative offset between the -0 angles of different detectors. 
The second is an overall offset in the focal plane reference 
frame and the frame on the sky. The latter is generated 
by any error in the calibration of the polarization angle of 
the instrument. 

Results from simulation of the if) systematics are shown 
in Figure[ni In the first two cases random 0.5 and 1 degree 
RMS -0 errors are added to each detector. This simulated 
fixed, random offsets in the relative polarization angles of 
the detectors. The results show the relative offsets con- 
tribute a comparable amount to the residuals as the nom- 
inal stepped mode case. 

In the remaining cases the same offset (0.25 degree, 0.5 
degree and 1 degree) is applied to all channels. This sim- 
ulates an overall calibration error in the half-wave plate 
angle. The results show this systematic gives a much 
larger contribution to the residual. For a 1 degree abso- 
lute offset the effect is as high as 20%. The 0.5 degree 
offsets is required having only a few percent residual ef- 
fect. Again, simulations consider only 8 pairs in a single 
column or 16 detectors total. The RMS result for the full 
focal plane should average down as ^/N , where N is the 
number of detectors. This factor has not been applied to 
the result. The calibration of the half- wave plate angles 
will be preformed pre-flight on the ground. The level of 
sub-percent precision required is not a difficult measure- 



ment and is made much easier by the compact optics and 
correspondingly close "far field" of Spider. 

Another systematic that is tested concerns pointing off- 
sets of crossed pairs. This occurs if the E-field distribution 
is not identical at the feed, or if there are polarization de- 
pendent properties in the optics. For these simulations one 
detector in a crossed pair (of the eight pairs) has a con- 
stant offset in RA and DEC of 1 or 4 arcminutes. Again, 
for the case of RA a factor of cos(DEC) is applied to the 
offset and hence represents the true offset on the sky. For 
Spider 40' beams, the corresponding A-B amplitudes are 
3.6% and 14% for 1 arcminute and 4 arcminute offsets. 
Residuals for beam offsets are plotted in Figure [TOl The 
effect is less than a few percent for £ < 80 for the largest 
offset case. 

4.4. Beam Systematics 

The Spider antenna array and optics define highly sym- 
metric beams on the sky. The beam pattern shown in 
Figure [2] is the feed beam pattern. The beam on the sky 
is influenced by the telescope. While the Spider telescope 
edge taper is modest, the beam on the sky will be more 
symmetric than the feed pattern shown here. In particu- 
lar, the visibly large and asymmetric lobes above will not 
propagate to the sky. The largest amplitude beam effect 
expected comes from reflections or "ghosting" in the Spi- 
der optics. Ghosting is common in refractive optics and 
results from unintended multiple reflections in the optics. 
The effect is a smaller amplitude beam image which is mir- 
rored with respect to the pixel position from the centre of 
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Fig. 10. — For these simulations one detector in a crossed-pair (of the eight pairs) has a constant offset in RA and DEC of 1 or 4 
arcminutes. For Spider 40' beams, the corresponding A-B amplitudes are 3.6% and 14% for 1 arcminute and 4 arcminute offsets. 
Residuals for beam offsets are plotted in Figure [TOl The effect is less than a few percent for ^ < 80 for the largest offset case. 
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Fig. 11. — The residuals from beam ghosting in the Spider optics. The ghosting is added in terms of a percent contamination to 
the nominal time-stream. Final maps are then reconstructed assuming no reflection contamination with no attempt to correct for 
the image distortion. A 10% contamination in the TOD yields a BB fractional residual as high as 20%. For 5% ghosting in the 
TOD this effect is already down by more than half and is negligible for 1% contamination in the TOD. 
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the focal plan. 

The ghosting is simulated by summing two time-streams 
from two beams. One of the time-streams is constructed 
using offset pointing from the ghost beam. The offset 
pointing of the ghost is determined in instrument coor- 
dinates (elevation and azimuth) , hence the ip angle for the 
ghost pointing is calculated appropriately as the final pro- 
jection of the orientation of a detector on the sky. The 
second time-stream is constructed from the pointing from 
the nominal beam. The two time-streams are summed 
with various weighting schemes depending on the ghost 
beam contaminatiorFi. 

The residuals from beam ghosting are summarized in 
Figure [11] For each case the ghosting effect is added in 
the time-stream in terms of a percent contamination added 
to the nominal time-stream. Final maps are then recon- 
structed assuming no reflection contamination; there is no 
attempt to correct for the image distortion. With 10% 
contamination in the TOD the effect in the BB residual is 
as much as 20%. For 5% ghosting in the TOD this effect is 
already down by more than half and is negligible for 1% in 
the TOD. Since simulations are done only for a single row 
(of 8 pairs of detectors) the distance between the original 
and ghost images is, on average, smaller for this row than 
for any other. It is therefore worth noting that the full 
focal plane may show a larger effect than that simulated 
here. 

4.5. Calibration Drift 

Diurnal variations in the detector sensitivity will occur 
due to altitude-induced changes in background loading. 
These sensitivity changes will be tracked using 4K semi- 
conductor emitters (fired intermittently) similar to those 
used in BOOMERANG flights. For BOOMERANG 
2003 flight, the gain drift for each individual detector was 
determined by the uncertainty on an individual calibration 
pulse and was measured to be 0.05%. With an improved 
version of the calibration lamp and Spider^ s higher sensi- 
tivity, the expected uncertainty is ~ 0.01%. Knowledge of 
the detector model (determined in pre-flight testing) will 
also allow the calculation of the sensitivity for any given 
operating point. 

Figure [T^ shows the effect of uncorrected calibration 
drift. Two cases are considered. For the first case the 
calibration drift is the same for all detectors, changing on 
a diurnal timescale with a maximum amplitude of 3%. For 
this case the resulting residual is small; less than a percent 
at all scales. 

For the second case the gain drift for all detectors are 
of the same amplitude and with a 24 hour period, but 
each of the 16 detectors have a gain drift with a different 
phase. Thus at any given time sample the calibration fac- 
tor within a pair of detectors will be different. Simulations 
consider only 8 pairs in a single column or 16 detectors to- 
tal. Again, the RMS result for the full focal plane should 
average down as -y/ZV, where N is the number of detec- 



tors. For this case a factor of v^/Vl024 (assuming 1024 
CMB science channels) has been applied to each residual 
result. The BB fractional residual is as much as 50% for a 
gain drift amplitude of 3%. For a 0.5% amplitude drift the 
residual drops below the 5% level. As with all simulations 
in this work there has been no attempt to correct for the 
gain systematic. For BOOMERANG 2003, the final re l- 
ative calibration uncertainty was 0.4% (iMasi et all 120061) . 
Spider is expected to achieve an uncertainty of 0.1% or less 
which will be more than adequate to meet science goals. 

5. CONCLUSIONS 

The results from Section 0] are summarized in terms of 
experimental specifications in Table [D While results are 
Spider-specific the order of magnitude of various effects 
can be translated to other CMB polarization experiments. 
The RMS B-mode signal for r = 0.1 is roughly 10000 nK^ 
(-1436 nK^ at ^ = 8 and -7379 nK^ at £ = 80), and 
scales linearly with r. Experimental specifications are set 
by limiting the allowed systematic residual level to a factor 
of — 10 smaller than the B-mode signal for r — 0.01. 

While the simulations were signal-only the impact of 
large low frequency detector noise (1/f noise) is reflected 
in the large scale degradation of the B-mode signal for the 
stepped half- wave plate mode of operation. Rapid, contin- 
uous half-wave plate modulation mitigates this effect en- 
tirely. Rapid, continuous gondola rotation also works but 
only with iterative map-making which accurately recovers 
the larger scale signal. 

It is important to note that the effects studied in Sec- 
tions 14.11 and Sections 14.21 (naive versus iterated maps, 
spinning more slowly, stepping half-wave plate versus 
spinning half-wave plate) will degrade the signal-to-noise 
achieved on the bandpowers. These effects differ from the 
systematics studied in Section 14.31 to Section 14.51 (point- 
ing reconstruction errors, polarization angle uncertainty, 
uncorrected ghosting, uncorrected gain drifts) which will 
ultimately bias the final result. The requirements on the 
biasing effects are more difficult treat than the signal-to- 
noise issues. 

The impact of systematics on B-mod e pola rimeter ex- 
periments is also discussed in iHu et al.] (|2003l) and more 
recently lO'Dea et al.l (|2007( ). where analytical methods 
are used for calculating the B-mode spectrum bias. The 
results are useful for setting experimental "benchmark pa- 
rameters" at the very earliest phases of instrument design. 
This work goes a step further by considering the impact 
of systematics in the map/time-domain; A necessary step 
in the evolution of an experiment which aims to measure 
the tiny primordial, gravity wave signal. 

This research used the McKenzie cluster at CITA, 
funded by the Canada Foundation for Innovation. Some 
of the results in thi s paper hav e been derived using the 
HEALPix pack age Gorski et a l. (2005)] as well as the 
FFTW package [Frigo fc Joh~n (2005, )]. 



° In order to completely model the polarization of the ghost, a half-wave plate angle dependency should be included in the ratio of ghost 
amplitude to main-beam response. This is not done here but is left for future work. 



14 



MacTavish et al. 



Calibration Drift 




& 



Fig. 12. — The effect of uncorrected calibration drift. In the first case the cahbration drift is the same for all detectors, changing 
on a diurnal timescale with a maximum amplitude of 3%. For the second case the gain drift for all detectors are of the same 
amplitude and diurnal but out of phase. For this case, the BB fractional residual is as much as 50% for a gain drift amplitude 
of 3% and drops below 5% for a 0.5% amplitude gain drift. Simulations consider only 8 pairs in a single column or 16 detectors 
total. The RMS result for the full focal pla ne wil l average down as y/N, where is the number of detectors. For the second case 
(out of phase gain drifts) a factor of \/T6/-\/1024 (assuming 1024 CMB science channels) has been applied to each residual result. 



Systematic 


Experimental Spec. 


Comments 


Receiver 




for llOdps 


1 /f knee 


< 200 mHz 


gondola spin 


Receiver 




for 36dps 


1/f knee 


< 100 mHz 


gondola spin 


Pointing 




sufficient for 


Jitter 


< 10' 


£<50 


Absolute 






Pol. Angle Offset 


< 0.25 deg. 




Relative 






Pol. Angle Offset 


< 1 deg. 




Knowledge of 




sufficient for 


Beam Centroids 


< 1' 


^ < 30 


Optical 
Ghosting 


< 2% 


% TOD 
contamination 


Calibration 






Drift 


< 3.0% 


in phase 


Calibration 






Drift 


< 0.1% 


out of phase 



Table 2 

Summary of experimental specifications based on simulation results. Realistic-amplitude, time-varying systematics are injected 
in the simulated time streams. Maps are reconstructed without any attempt to correct for the systematic errors. Experimental 

specifications are set by limiting the allowed systematic residual level to a factor of ~ 10 smaller than the B-mode signal for 
r = 0.01. The nominal operating mode is a 36 dps gondola spin rate, with the half-wave plate stepping 22.5° once per day, with 
10 iterations (sufficient to recover the residual levels of the cont inuously-rotating half-wave plate case) of the map-maker, a 

Jacobi iterative solver (jJones et al.ll2007l ). 
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